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Basic ideas and results which characterize quantum diffusion of defects in quantum crystals like 
solid helium as a new phenomenon are presented. Quantum effects in such media lead to a der- 
ealization of point defects (vacancies, impurities etc.) and they turn into quasiparticles of a new 
type — defectons, which are characterized not by their position in the crystal lattice but by their 
quasimomentum and dispersion law. Defecton-defecton and defecton-phonon scattering are consid- 
ered and an interpolation formula for the diffusion coefficient valid in all interesting temperature 
and concentration regions is presented. A comparison with the experimental data is made. Some 
alternative points of view are discussed in detail and the inconsistency of the Kisvarsanyi-Sullivan 
theory is shown. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The basic characteristic of a quantum crystal is the 
presence of a large amplitude of atomic zero-point vibra- 
tions. This means that the overlap of the wave functions 
of the nearest neighbours is not negligible and leads to 
considerable tunnelling transitions to the nearest lattice 
sites. If a point defect (vacancy, impurity, interstitial) 
appears in such a system, the crystal state becomes de- 
generate with respect to the defect position due to the 
identity of lattice sites and the defect delocalizes and 
turns into a quasiparticle - defecton. If the wave func- 
tion of the crystal with a localized defect in the lattice 
site r is iprix), then the right wave function of the crystal, 
according to the Bloch theorem, is Q 



^Vr(a;)exp{i(kr- e(k)i/fi)} (1) 



where x is a 37V component vector which consists of all 
atomic coordinates, and k is the defecton quasiwavevec- 
tor. The spectrum is given by the defecton dispersion law 
e(k). In a tight binding approximation it has the form 



e(k) = $ + Aw(k) 



(2) 



where $ is the activation energy of a localized defect, A is 
the tunnelling probability amplitude (exchange integral), 
and w(k) is a dimensionless function which depends on 
the symmetry of the crystal lattice. For a simple cubic 
lattice the dispersion law is of the form 



e(k) = $ + v4(cos k x a + cos k y a + cos k z a) 

and for body-centered cubic lattice (b.c.c.) 
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where a is the lattice constant. For defects in solid he- 
lium, the quantity $ is of the order of several degrees. In 
case of h.c.p. lattice, there are two atoms per unit cell 



and the dispersion law consists of two branches called 
conditionally acoustic and optical respectively [HSU 
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An important characteristic of a defecton is its energy 
band width A. For simple lattices A = zA, where z is 
the number of nearest neighbours in the lattice. For the 
h.c.p. lattice, the dispersion law (J5J is written in a form 
which fixes the bottom of the acoustic band at k = 0, so 
that e(0) = en is the lowest band level. The band width of 
the acoustic branch is A ac = 7 A while that of the optical 
one it is A opt = 5 A. The bottom of the optical branch lies 
higher by 3A A typical value of A for 3 He impurities in 
4 He is A - IQ-^K, while for vacancions A is of the order 

1 de 

of several degrees. The defecton velocity is v = — — 

h OK 

and may be estimated by the relation v ~ oA/fi. The 
effective mass for simple lattices is m* = ti 2 /Aa 2 while 
for the h.c.p. dispersion law (JSJ) one has m1 cp = h 2 /2Aa 2 . 
Another important quantity is the derealization time To 
which may be estimated by the relation 0, EJ 



T Q ~ A/h. 



(7) 



Note, that in this relation the characteristic energy is the 
exchange integral A, not the band width A. 

Before going on with defecton diffusion let us note that 
there are crystals in which zero-point vibrations are large 
only for certain degrees of freedom and hence the quan- 
tum properties manifest themselves with respect to them. 
The class of quantum crystals includes therefore heavy 
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metals (palladium, niobium etc.) in which the behaviour 
of light hydrogen impurities (protons) is quantum, solid 
hydrogen with impurities of hydrogen atoms, metallic hy- 
drogen, etc. 

From a historical point of view, however, physics of 
quantum crystals is above all bound up with solid he- 
lium. Helium has been chosen for a basic system when 
examining the ideas and results of the theory owing to 
a number of circumstances. Solid helium is the best 
investigated quantum crystal. Helium crystals are ex- 
tremely pure - all other elements solidify at higher tem- 
peratures, so the only impurities are isotopes and can be 
well controlled in wide concentration ranges. The two 
helium isotopes obey two quantum statistics - those of 
Bose-Einstein and Fermi-Dirac. In addition, each iso- 
tope crystallizes in three different modifications. In such 
a way solid helium offers six quantum crystals, composed 
of two isotopes only. Phase diagrams for helium isotopes 
can be found, e.g. in 

Any correct theoretical model has to include not only 
a description of the structure and building particles, but 
also the typical kinds of motion in it, or, in other words, 
the elementary excitations as elementary carriers of the 
motions. In this sense crystals of solid 4 He with small 
concentrations of isotopic impurities and vacancies may 
serve as examples of really existing model systems, in 
which the only elementary excitations are phonons and 
defectons. In addition, the low value of the melting tem- 
perature compared to the Debye temperature leads to 
a simple linear dispersion law for phonons in the whole 
range of existence of the solid phase. It is worth noting 
that the new (tunnel) kind of motion presented by the 
defecton quasiparticles is the reason to call the state of 
quantum crystals the a new state of the matter. 

For 3 He atoms which have spin equal to 1/2 and rel- 
atively large magnetic moment an exchange of atoms is 
equivalent to a spin exchange (the interatomic potential 
practically does not depend on spin) and therefore may 
be observed experimentally, e.g. by means of nuclear 
magnetic resonance (NMR) 9]. 

The presence of two types of crystals composed of 
Fermi- and Bose- atoms puts forward the problem of 
defecton statistics. It can be shown that in pure crys- 
tals with small defect concentration vacancions obey the 
statistics of host atoms. Impurities of 3 He in crystals 
of 4 He as well as impurities of 4 He in crystals of 3 He 
turn into Fermi-quasiparticles. In case of mixed crystals 
a vacancy has to exchange places with different atoms 
( 3 He or 4 He) that obey different statistics, therefore, the 
vacancion statistics can differ from the known statistics 
of Bose-Einstein or Fermi-Dirac. More information on 
defecton statistics can be found elsewhere 0, ^} . The 
degeneracy temperature T is extremely small: 

T ~ Ax 2/3 (8) 

where x is the fractional concentration of defects. For 
vacancions with a typical concentration x ~ 10~ 3 the 
degeneracy temperature is Tq ~ 10~ 2 — 10~ 3 K. That 



is why in many cases defectons may be considered as a 
Boltzmann gas. It is important to note, however, that the 
role of statistics increases with decreasing concentration. 
Indeed, the interaction energy Ei nt between point defects 
in solids is inversely proportional to the third power of the 
spacing r and hence is proportional to the concentration 
x: 

Eint = ^=V x (9) 

where Vq is the elastic interaction constant which typical 
value is of the order of 10~ 2 K . Taking into account JS}, 
one obtains 



Therefore, at low enough concentrations the interaction 
energy E turns out to be negligible compared to the 
"quantum statistics energy" To. The critical concentra- 
tion for a gas of vacancies is ~ 1CP 1 and is much greater 
than the usual vacancy concentration in the experiments. 

II. HISTORICAL REVIEW 

A phenomenological theory of defects in quantum crys- 
tals was proposed first by A.F.Andreev and I.M.Lifshitz 
1 1 1 1] in 1969. The qualitative estimations concerning dif- 
fusion in that work differ from the quantitative ones more 
than 7 orders of magnitude and are only of historical 
interest now. Simultaneously, the author of this work 
reported at the Second Soviet Conference on Theory of 
Solid State in Moscow his microscopic theory of vacan- 
cions and impuritons (published as a paper in 1970) 
where the defecton parameters were estimated and the 
temperature dependence of the defecton mean free path 
was obtained. 

At the same time Guyer and Zane ^2 published their 
perturbation theory of the spectrum of impuritons (called 
mass-fluctuation waves) and showed the possibility of ex- 
perimental observation by the method of nuclear mag- 
netic resonance. 

It is worth noting that there was a significant inter- 
est in quantum crystals at those years but there were no 
experimental evidences on quantum diffusion. Moreover, 
some results on the diffusion of point defects [l3L Iblj 
were considered as arguments in favour of the classical 
activation mechanism. There was not any clear concep- 
tion of the correlation between the quantum and clas- 
sical weights in the phenomena considered. This pre- 
vented often even the correct formulation of the prob- 
lems. Instructive in this aspect is the work of C.P.Flynn 
and A.M.Stoneham [15J on the diffusion of light inter- 
stitials. The authors have practically neglected the co- 
herent quantum diffusion stressing the role of phonons 
and deformation around the defect. This was a re- 
sult of their conception that the incoherent (by means 
of phonons) transition corresponded to a macroscopic 
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number of channels, while the coherent tunnelling cor- 
responded to one. In this way, an implausible estimation 
has been made for the temperature region in which coher- 
ent diffusion could be observed, given by the inequality 
T < 10- 17 K. 

The assumption that defects turn into quasiparticles 
means that (at some conditions considered later on) the 
gas kinetic model is applicable and hence the diffusion 
coefficient has to be inversely proportional to the con- 
centration, i.e. 



D ~ Aa 2 /Hxa 



(11) 



where a is the defecton-defecton cross-section (in units 
a 2 ). The temperature dependence predicted was 
of the form 



Dn 



A 



*-5- < 12) 



The ninth power of temperature is easy to understand 
(for details see Refs. [lj, fl8^> . The motion of a defect is 
determined by the number of phonons (N p h ~ T 3 at 
low temperatures) and, by the scattering of long wave 
phonons on defects (Rayleigh scattering) with a cross- 
section a ~ fc 4 ~ T 4 ; finally, the motion of a "heavy and 
slow" defecton in a gas of "light" phonons depends on the 
second power of the phonon wave vector (~ k 2 ~ T 2 ). 
Hence, D ~ l/(N ph aT 2 ) ~ T" 9 . It is worth noting 
that the role of the characteristic temperature is played 
by 9p, not by the Debye temperature Qu- These two 

temperatures differ significantly: 9 p = — rs — — and 

2a 8 

lead to a numerical difference of the order of 8 7 ~ 10 6 
times. 

It could be expected that such a rapid increase of 
the diffusion with decreasing temperature can easily be 
observed experimentally. However, the experiments on 
both vacancy |13| and impurity diffusion carried out 
in He with 0.75 -r- 2% impurity concentration indicated 
typical classical exponential decreasing down to the phase 
separation temperature. 

The first evidence of a quasiparticle behaviour was 
reported by R.G.Richards, J.Pope and A.Widom pcj . 
They confirmed the concentration dependence ~ 1/x. 
However for the tunnelling frequency A they obtained 
the value A ~ 10~ 7 K, which is three orders of magni- 
tude smaller than the exchange integral J33 in pure 3 He. 
The interpretation of such a value ran into difficulties be- 
cause there were no physical reasons for such dramatic 
difference in the exchange integrals. It became clear that 
additional information is necessary and that it could be 
obtained from the temperature dependence (|12l) . An in- 
creasing diffusion with decreasing temperature was first 
observed by the Kharkov group [2l|,|22] in 1973. Surpris- 
ingly, the treatment of the experimental data using the 
theory of Andreev and Lifshitz ^lj gave for the exchange 
integral A a value of the same order of magnitude as that 
obtained in 1201. 



The plain model of free quasiparticles became ques- 
tionable and in 1974 Landesman and Winter (23 and 
later on Huang et al. [24[ proposed an alternative model 
leading to another concentration dependence: 



D 



Ja 2 



K 



(13) 



where J' is the exchange integral in the pure crystal, 
and K is a phenomenological parameter describing the 
repulsion of two atoms in one lattice site. 

The authors argued that this dependence is valid for 
concentrations x > 10~ 6 , i.e. in the whole experimen- 
tally accessible region. Unfortunately, the experimental 
precision was not enough to distinguish the small differ- 
ence between the concentration dependencies in (|llfl and 
(|13fl . especially taking into account that both A and J 
were unknown. 

In order to solve the paradox the author treated the 
experimental data of the Kharkov group using his theory 
[l7| . As a result, the exchange integral A was found to 
be only slightly smaller than J 33 , the band width was 
estimated to be A ~ 10 K and the prediction of the gi- 
ant defecton-defecton cross-sections in pj was confirmed. 
The existence of the new quasiparticles (defectons) be- 
came an immutable fact and the discovery of a new phe- 
nomenon called quantum diffusion was acknowledged. 



III. DEFECTON-DEFECTON SCATTERING 

The only scattering processes which determine the 
defecton diffusion are collisions with defectons and 
phonons. If the temperature is small enough (the exact 
conditions will be derived below) the number of phonons 
is negligible, and the defecton free path is controlled by 
defecton-defecton interaction. 

Defecton-defecton scattering possesses a series ofpecu- 
liarities relevant to the energy band width P, |3, E3) 13 • 
First of all, the quasiparticle energy cannot vary during 
the motion more than the energy band width A. This 
means that the quasiparticle is not able to overcome bar- 
riers higher than A. In case of impuritons in helium, a 
typical value of the deformation energy caused by a va- 
cancy or impurity is of the order of 10~ 2 K and is much 
larger than Ai mp ~ 10 -4 K. As a result, the quasiparti- 
cle cannot reach the core of the vacancy (or interstitial) 
and exchange place with it (the later case of vacancy- 
stimulated diffusion has been discussed in Sec. lVII|) . This 
is important because the exact form of the deformation 
in the core is complicated and not known, while the be- 
haviour of the deformation field at large distances is com- 
mon for wide classes of solid media. The later allows to 
study the defecton scattering in a quite general form. 

Point defects in crystals interact by the law 



tp(r) = V T(n) (- 
V r 



(14) 
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where r(n) is a function of the direction n, which van- 
ishes after angle averaging and is of the order of unity. 
In an isotropic media 



9 (r) = V (- 
V r 



(15) 



Let us consider first the features of the defecton scatter- 
ing in a potential of a more general form: 



fix) = Vo (- 
V r 



(16) 



The small value of the impuriton energy band width leads 
to a peculiar confinement in areas with strong field gra- 
dients. This effect takes place both in attractive and in 
repulsive potentials. Its physical reason is that the tran- 
sition of the quasiparticle into a neighbouring lattice site 
cannot take more energy than the band width. There- 
fore, the condition for such a confinement may be written 
in the form 



Vtp | a > A. 



(17) 



Substituting here tp from (|16fl and replacing the mean 
distance r between defects by their concentration x ~ 
one obtains that all defects has to be localized 



a/r l 



111 at concentrations x > x r , where 



3 

1 + 1! 



A 

nV 



In the case of isotopic impurities in helium (n 
quantity x c is 



A 
3^ 



1 - 2%. 



(18) 



3) the 



(19) 



Due to the small value of the defecton band width large 
regions around the impurities are unscalable for quasi- 
particles [lLITs|. Their linear size Ro may be estimated 
from the equation 



p(iZo) = A, Ro = a 



l/r, 



(20) 



Hence, the effective cross-section a for impuritons in 
the potential Ijltifl is 



Vo 



2/r, 



(21) 



These estimations may be improved taking into ac- 
count, that the experiments are performed at temper- 
ature T ^> A and, therefore, the conditions for the 
quasi-classical approximation are satisfied: 



(feet 



,71-2 



A _ A 
Vb ~ ~zVo" 



(22) 



In this case the cross-section is given by the expression 



a « 2na 



Vo_ 
Aka 



V 



(23) 



(States with ka > 1 are activated due to the high tem- 
perature) . 

In the most interesting case, when n — 3, the cross- 
section 



,v 

A 



(24) 



can reach giant values a ~ 10 2 a 2 . Eq. Il2,'jll shows that the 
effective cross-section in the gas approximation (nonin- 
teracting quasiparticles) considered is inversely propor- 
tional to the band width A. The same dependence on 
A follows also from the theories of Landesman [25| and 
Yamashita .26] for interacting quasiparticles. This co- 
incidence is evidently the explanation of the fact that 
the dependence of the diffusion coefficient on the molar 
volume ( "the Gruneisen parameter" ) does not vary with 
decreasing concentration [22| (see below). 

The condition 122|) can be satisfied at lower tempera- 
tures as well if the following inequalities take place 



A(A/V )— «T« A 



(25) 



Then the quasipartitlc thermal energy may be evaluated 
from the relation h 2 k 2 /2m* w 3T/2. With to* ps H 2 /Aa 2 
this yields ka ~ (T/A) 1 / 2 , and the cross-section (|2l5|) can 
be written in the form 



V 2 
TA 



Vo 



71=3 



(26) 



It increases with decreasing temperature up to the value 



l {V0/Ay 



KVo/Af 



(27) 



The opposite situation which corresponds to the scatter- 
ing of a slow particle by a rapidly decreasing potential 
barrier [n = 6) can take place for vacancions at T <C A. 
Then the cross-section does not depend on velocity, and 
hence, is independent of temperature: 



<t ^ 3.657ra 2 ( — %■ 
2A 



(28) 



This expression fits well to that of Eq. (|27|l at n = 6. 
If vacancion-vacancion scattering is under consideration, 
then Vo ~ ms 2 , A = IK, and hence, a = 50a 2 . 

In the consideration above we did not take into ac- 
count the role of the anisotropy given by the factor r(n) 
in the potential (2J. It is obvious that such a sign- 
variable function with zero averaged value can only re- 
duce the cross-section. According to Slusarev et al. 
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the anisotropy reduces the total cross-section by a fac- 
tor 0,105 (cubic lattices), and 0,208 (h.c.p. crystals). 
The reduction for the transport cross-section is about 4-5 
times. However, the exchange integral A as well the de- 
fecton band width A evaluated in the same work turned 
out to be more than an order of magnitude less compared 
to the results adopted now. 

Situations where A ^> Vq are realized in vacancion 
scattering by isotopic impurities. In such cases Born ap- 
proximation may be used and the problem can fully be 
solved in an analytical form with an account of the ex- 
act form of the function r(n) 0,0]. It can be shown 
that the scattering amplitude F(n) has the same angular 
dependence as the interaction potential: F(n) = Ar(n). 
For a simple cubic lattice A = (2a / 15) (Vo / A) , and there- 
fore, 



da = —a 2 (V /A) 2 T 2 (n)dn, 
225 



(29) 



where r(n) = 



3/5. 



One has to keep in mind, however, that the interaction 
potential has the form (| 1 41) only at large distances and 
the effect of the scattering by the core needs for wide- 
band vacancions a special consideration. 

Consider briefly the role of the defecton statistics. If 
defectons obey Fermi-Dirac statistics, then their quasi- 
momenta has to be near to the value of the Fermi- 
quasimomentum both before and after collision. As far as 
the cross-section is proportional to the final states num- 
ber, it is proportional to temperature T. Therefore, 



op = a'T/ep. 



The quantity a' for the potential with n — 3 in 
quasiclassical approximation is of the form 0, l29| 



a' = 2tt 



2 Voa 

hvp 



(31) 



where vp is the Fermi-velocity. Therefore, the cross- 
section depends not only on the ratio Vo/A, but also on 
the concentration: 



<tf~T^-, T<Ax 2 '\ 
A A x 



(32) 



The condition for the applicability of the expressions 131|) 
and (|32(l follows from (1221) where one must take into ac- 
count, that kpa ~ x 1 / 3 : 



x > (A/V ) 



(33) 



It is worth noting also that the ratio Vo/A is strongly 
dependent on pressure, and increases rapidly with de- 
creasing molar volume. 



IV. LOW TEMPERATURE DIFFUSION 

Let us now turn to the diffusion. As follows from the 
jas kinetics theory the diffusion coefficient is equal to 

1, 



D 



(34) 



where the mean free path is 

I = a/V2xa 



(35) 



This approximation is valid if the mean distance between 
quasiparticles, £ -1 / 3 , is much larger than the "radius" of 
-^/er, i.e. if x^ 1 ^ 3 ^> \fa or (what 



the cross-section R 
is the same): 



x « 1/a 3 / 2 . 



(36) 



The quasiparticle mean velocity v depends, in general, 
on the temperature. If T 3> A which is the case realized 
in experiments on impurity diffusion in helium, then the 
levels till the top of the energy band are populated. Sacco 
and Widom [3(| calculated the band-averaged value of 
the velocity squared in a h.c.p. lattice and obtained (in 
our notation) 



1 . 2 . si/2 3 a Q A 
vsw = {<v >) * = 



a/2 h 



3 A ~ A 

27I S ^ V (37) 



In principle, this value can be used to estimate the order 
of magnitude of the diffusion coefficient l|34|l . However, 
it is well known (see e.g. Ref. |3l|) that the correct value 
of the velocity is closer to its highest value which can be 
evaluated from the relation m*v 2 /2 = A: 



1 A 

/2z0 P 



2,45s- 



A 



(z = l2). 



(38) 



(30) The diffusion coefficient is therefore (comp. to (lll(l ): 



Dn = 



1 as A 
xa 8 p 



1 as A 
V3xa 6 P 



(z = 12). 



Substituting here a from yields 



D 



A 2 



A 2 



677 y/z x VoOp 3iryfz hx Vo 



(39) 



(40) 



In order for the defectons to be good quasiparticles, their 
mean free path should be larger than the lattice constant, 
i.e. the coherency condition xa < 1 should be satisfied. 
This leads to the following restriction on the defecton 
concentration: 



x < xo 



(41) 



As we showed in the previous section the cross-sections 
a can be extremely large. That is why the condition (|41() 
is weaker than (|36|l . Hence, the gas-kinetics formula (|39|) 
is valid for concentrations 



x < a 



- 3 / 2 < a- 1 = xo. 



(42) 



The concentration xo turns out to be smaller than the 
critical concentration x c 1)19(1. This means that the lo- 
calization process may be considered in terms of quasi- 
particles. Making use of the results obtained we may 
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distinguish the following concentration regions depend- 
ing on the powers of the small parameter £ = A/Vq : 



I. x < £ 3/,l_1 < 10" 

II. C 37 "- 1 < x < (£/n) 3 /" +1 
quasiparticles 



gas approximation 



10- 



interacting 



localized defectons 



III. {i/nf/ n+1 < x 
(The numerical values are calculated for n = 3) 

Let us consider now the diffusion in the second region. 
In this case the mean free path / is controlled by the scat- 
tering on the defecton band boundaries (Fig. ??.), i.e. by 
the condition I \ \j(p |w A. This means that the vari- 
ation of the defecton energy along a distance I should 
not exceed the band width A. The corresponding dif- 
fusion coefficient in this region is therefore of the form 



D 



A 2 



V oX (n+l)/3 



(43) 



It is worth noting that the cross-sections both in re- 
gion II, and in region I depend on A. This dependence 
is important not only for choosing suitable defect con- 
centration and perfection of the crystal samples. The 
knowledge of the diffusion coefficient as a function of A 
allows to find the dependence of D on the molar volume 
V m . Since this dependence is strong and different in dif- 
ferent concentration regions, one may use it to detect a 
change of the scattering mechanism. 

At temperatures lower than the degeneracy tempera- 
ture To, diffusion depends on defecton statistics [23. If 
the quasiparticles obey Fermi-statistics, then the gas ap- 
proximation yields 

Dp = -l F v F , where v F = (37T 2 ) 1 / 3 — x 1/3 (44) 
3 n 

and the mean free path is 



Of course, the increase of l F and Dp with decreasing 
temperature is limited by the finite size of the sample, or 
by scattering on other crystal defects. 

In the case of Bose-defectons the diffusion coefficient 
may be estimated taking into account that the "actual" 
concentration depends only on temperature: 



1 (T 



* W 6 I 



3/2 



and hence, the mean free path is 



6 a 



V2a \TJ 



The mean velocity is 

uwl,75j( AT) 1 / 2 , 
h 

The diffusion coefficient is therefore 
D«1,26. M * 



a 6 P T 



If the mean free path turns out to be larger than the 
sample size L, then 



D 



^{AT) 1 ' 2 . 



One has to keep in mind, however, that those effects 
can be violated by a superfiuid current, which can take 
place even at small chemical potential gradient. 

Apparently, the diffusion at T < Tq is not observable 
in helium due to the rapid decreasing of vacancion con- 
centration (except of zero-point vacancions, not found 
yet) and the phase separation in the case of impuritons. 
To observe it one may have to look for another quantum 
system. 



\f2ax e ff 

Since only defectons near the Fermi-surface take part, 
the "active concentration" is equal to 



x eff = x- 



T 
e F 



and the effective cross-section is given by Eqs. i[30 |l — (|32 |l . 
As a result, the mean free path and the diffusion coef- 
ficient are multiplied by an additional factor (e F /T) 2 . 
This yields 



D F = 



1 av F (e F \ 2 

3 X(7q 



(45) 



Taking into account (|31|) one may rewrite Eq. 145|l in 
the form 



37T 2 A 4 



(46) 



V. PHONON CONTROLLED DIFFUSION 

At absolute zero temperature a defecton moves 
through a crystal like a free particle without dissipation. 
In view of its low velocity compared with the sound one, 
the gas of zero-point phonons has time to adjust itself to 
the defecton which moves together with an adiabatically 
adapted phonon "cloud" that strongly influences only its 
effective mass. Note, that null-phonons do not change 
the periodicity of the lattice. 

With rising temperature, scattering of thermal 
phonons by the defecton sets in. Since the defecton mass 
is large and its lifetime is much longer that the relaxation 
time of the phonon gas the cross section for scattering by 
a defecton can be replaced in zeroth approximation with 
respect to the small parameter (v/s) 2 by the cross sec- 
tion for scattering by an immobile defect. As it has been 
shown by I.M.Lifshitz [3^ (see also 0111]) such a prob- 
lem can be solved in an analytical form pi 0, Q • 
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The differential cross-section for the scattering of a 
phonon of the j-th mode with polarization v as a result 
of which a phonon of the i-th mode with a polarization 
v' is produced is given by the expression 0: 



aa jW,ik' 



( v 



V 47TTOS,; 



(<?rr 



£kRA^,k|R'> ,/o'M7) 



RR> 



where V is the lattice cell volume, s, is the sound ve- 
locity of the corresponding phonon branch, is the po- 
larization vector and dfl' — 2n sin ada, a being the an- 
gle between the directions of the incident and scattered 
phonons. The matrix A RR , (r) describes the change of 
the elasticity matrix caused by the defect (situated in 
the lattice site r). Averaging over the modes and polar- 
izations of the incident phonons and summing over the 
modes and polarizations of the scattered ones yield 



d<T k k> = ( — ) fc 4 S nn /dr2' 



(48) 



where n = k/fe, n' = k'/A; 

V 



12m 2 s 4 



W^nRA^n'R'j 

fiv IrR' J 



Note, that the averaged value of the sound velocity 
s = (s ; ~ 4 + 2sJ 4 )~ 1 / 4 does not coincide with the Debye 
velocity, even if does not differ significantly. 

Any further application of the formula (|48|) requires 
some model assumptions regarding the matrix compo- 
nents AJ^, . They have to obey, however, some relations 
that follow from the behaviour of the forces in a crystal 
lattice 36] . In the nearest neighbour approximation this 
allows to express all matrix elements by one of them (e.g. 
Aoo) 111- In case of a vacancy in a simple cubic lattice 
Aoo = — 6ms 2 /a 2 and no free parameters are required. 
In this case V = a 3 , T,„„ 



da 



kk' 



2tt 



(nn') 2 and 
(fca) 4 (nn') 2 c«y. 



(49) 



In case of impurity, this quantity should be multiplied by 
a factor a = (A ™ pu "* y / 1 ^ ) cancv ) 2 Note, that during all 
the consideration we have not supposed the perturbation 
caused by the defect small. This is the advantage of the 
Lifshitz method used. Let us call attention to the fol- 
lowing peculiarities when considering isotope defects. In 
a classical crystal isotope defects interact with the host 
atoms in approximately the same way and the main dif- 
ference consists in their mass. This leads to a correction 
factor cr = (1 - m host /m isotope ) 2 . For 3 He impurities 
in solid 4 He ao — 1/9. In case of quantum crystals iso- 
topes differ not only in mass but also in their zero-point 



vibrations affecting in this way the local values of the 
elasticity module 0. An attempt to estimate this effect 
was made by Slusarev et al. |37|. They argued that this 
effect could be larger than that of the mass-difference 
and that we had not taken it into account. The later 
statement is not acceptable as seen by the consideration 
above. In our short letter the principal value of the 
cross-section (given by the mass-difference) was used in 
order to avoid model considerations. We shall not dis- 
cuss here the model assumptions in [3?} in detail. How- 
ever, the correction proposed is given by the quantity 
d « 83/8- 1 « 0.23 in the expression (cf. Eq.(15) in 

Ref.|aa) 



Gi/4 



a 2 



1 



0.11 + 0.0145 



(50) 



with a = ao = 1/3. Obviously, the correction 0.0145 
could not be considered as 'larger' than 0.11. In our 
opinion next terms in Eq.(14) in [37j of the order of 
Gi/2 12 ~ 2 x 10~ 4 Gi are far from the model accuracy 
used. 

The diffusion coefficient may be calculated solving the 
Boltzmann equation, which is of the type of Fokker-Plank 
equation However, the specific features of the de- 

fecton motion in the phonon gas allows to solve the prob- 
lem without resorting to the kinetic equation. We follow 
here the method proposed in 0, 0] . 

Since the defecton velocity is much smaller than the 
sound one the main contribution is due to two-phonon 
processes. The energy and quasimomentum conservation 
laws for an individual collision are: 



e(p) + hsk = e(p') + Hsk' 



p + k = p' + k' + 2?rb 



(51) 



(52) 



where k and k' are the quasiwave vectors of the incident 
and scattered phonons respectively, e(p) - the defecton 
dispersion law, and b is an arbitrary vector of the recip- 
rocal lattice. The linear dispersion law for the phonons 
is used since the temperature is much smaller than the 
Debye one. The phonon gas can be considered as in equi- 
librium, and hence the mean value of the phonon quasi- 
momentum is hk cx T/s. Only phonons with quasiwave 
vectors near the centre of the Brillouin zone take part 
in collisions, and Umklapp processes may be neglected. 
The temperature, in general, can be higher or lower than 
the defecton band width A. In both cases the relative 
variation of the defecton energy and quasimomentum in 
a collision are small. If A <C T <C 9 P then Eqs. I|51l) and 
yield 

de 

Se = —Sp = ftv(p - p') = Sv(k-k') < 2vhk cx 2vT/s 



and hence, 



Se/A < T/0 p < 1 



(53) 
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In the same way 

| p - p' | | k - k' | 2T T 



< 



P 



P 



hsp paOp 



< 1, ti a aA/S 
(54) 

In the opposite case, when T < A < 9 p , the defecton 
dispersion law may be replaced by a quadratic one and 
the effective mass approximation can be used: 



h 2 p 2 

2m* 



Then the change of the phonon wave vector is 

k k 2 

X = k' - k = — p(n'-n) (1 - nn') (55) 

Pq Po 

where po = (m*s)/H = 28 p /(Aa) and n = k/fc, n' = k'/fc 
— the directions of the incident and scattered phonons. 
The wave vector transferred is 



q = fc(n - n') + x n ' ■ 



(56) 



The relative variations of energy and quasimomentum 
are, therefore, small in this case too 



6E k (AT) 1 / 2 

-— oc - oc — 

e(p) p Op 



q k 
< 1, - oc - < 1. 
P P 



As for the phonon quasimomentum, its direction changes 
greatly in each collision, but its magnitude remains prac- 
tically constant. The situation thus recalls the motion of 
a heavy particle in a gas of light particles. The difference 
lies in the dispersion laws and in the different statistics 
obeyed by the defectons and phonons. It is important 
also that the width of the defecton energy band is rela- 
tively small. 

Let us find first the diffusion coefficient in momentum 
space. Obviously, the change of the square of the wave 
vector per unit time is given by the integral 



<9 2 > = T^TT / d 3 kd<rwn(k)(k-k' 
(27r) d a 6 J 



(57) 



where n(k) — {exp(fofc/T) — 1} is the phonon distri- 
bution function, and dckk' is the differential cross sec- 
tion l]48p. At low temperatures (T <C 9 P ) the integral 
in Eq. I|57|) may be expressed in terms of the Riemann 
(■-function: 



a s ( T 



[Q ) = ot^r — 



7T 5 C(9) 

240C(8) 



1,27. 



(58) 



To calculate the diffusion coefficient in coordinate 
space we proceed as follows. We define the free path 
time of the defecton as a time during which the trans- 
ferred momentum squared becomes of the order of the 
square of the initial particle momentum: 



l_e(p) 
aero 9 p 



T 



T = h/A 



(59) 



Then the mean free path is the path traversed during 
the time r : 



/ = 



V2 a e(pf /2 (0^ 9 



oi <Jo Q V \J~A V T 



(60) 



and the diffusion coefficient in the phonon gas is 

*=^(?) 2 (D 9 <«> 

The presence of the ninth power the temperature has a 
lucid physical meaning - three degrees are connected with 
the number of phonons, four with the scattering cross sec- 
tion, and two with the "ineffectiveness" of the collisions 
(when a heavy particle is scattered by light ones). 

It is instructive to note once again that in all the equa- 
tions the characteristic temperature parameter is not the 
Debye temperature 6 as assumed in 0, 0, 0, but 
p sa 0/8. The effect of this is that the kinetic character- 
istics evaluated in the works mentioned differ from our 
results by a great facto r of the order of 10 8 . 

The mean value e(p) depends both on the temperature 
region and on t he q uasiparticle statistics. For Fermi- 
Dirac statistics, e(p) = cf oc Ax 2 / 3 , and 



D F oc -x 4 / 3 



2 /f) \ 9 

Vp 

T 



(TQ yup 

In the case of Bose-Einstcin statistics: 



(62) 



D B cx- ^ 



CO 



T 
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(63) 



Since the degeneracy temperature To oc Ax 2 / 3 is ex- 
tremely small, the inequality To <C T is practically satis- 
fied in the whole experimentally accessible region. If at 
the same time T < A, then Boltzmann statistics applies 
and the relation e(p) — 3T/2 may be used to estimate the 
averages over the temperature. In this case we obtain for 
the free path time and for the diffusion coefficient, re- 
spectively, 



7*«1,182L £ 

CO \T 



Dr 



0.59^| . (64) 



From the last expression and the Einstein relation we 
easily obtain the mobility b of the defecton in the phonon 

gas 



Db 
T 



1.18 



ha 



(65) 



Naturally, the same value of the mobility can be obtained 
by direct calculation of the force experienced by a quasi- 
particle moving in the gas. Indeed, in the system in which 
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the defecton is at rest, the phonon distribution function 
takes the form n(k) = n(e — fikv), and the force act- 
ing on the defecton can be expressed as the change of 
momentum per unit time: 



3sh 

(27T) 5 



J d 3 kda kk >n(k)k(n - n') = b x v. 



In first order in the small ratio v/s <C 1, the mobility 
calculated in this manner coincides with that of Eq. (|65l) . 

If T > A, then the band is filled uniformly and e(p) = 
A. In this region which is the most interesting one from 
the experimental point of view we have 



tq A 
aa P 



D i 



A 



(66) 



It follows from the deduction made that the results 
obtained are valid as far as the mean free path I is much 
longer than the lattice constant (I ^> a). This condition 
leads to the inequality T <Tk, where 



Th — On 



V2z A 



1/9 



(67) 



At temperatures T > T k the mean free path becomes 
smaller than the interatomic distances, and the free path 
time r gets shorter than the derealization time J7|). As 
a result, the defecton spends a greater part of the time 
within the cell, and only rarely does it execute individ- 
ual transitions (tunnel or activation) to a neighbouring 
equivalent position — the so called phonon stimulated 
localization takes place. At T — T k the diffusion coeffi- 
cient is 



D cx 



a 2 A 



(68) 



Its behaviour at higher temperatures calls for a special 
analysis. 

It has been reported in [3^ that in a case of light in- 
terstitials the law D ~ T~ 9 is valid up to the Debye 
temperature, independently of the small mean free path. 
Unfortunately, the melting point of solid helium is much 
lower the Debye temperature, and hence this prediction 
cannot be tested experimentally. On the other hand the 
consideration made in ^ij^ does not take into account 
that the activation mechanism may become essential far 
below the Debye temperature 0. 



VI. COMPARISON WITH EXPERIMENT. 
DETERMINATION OF THE IMPURITON 
CHARACTERISTICS. 

Let us summarize the results obtained and compare 
them with the available experimental data on diffusion 
of He impurities in solid He. 



The most important processes which determine quan- 
tum diffusion in helium are: i) scattering by defects, 
leading to a temperature independent diffusion, and 
ii) defecton-phonon scattering which determine the tem- 
perature dependence. In the most interesting case when 
T ^> A the diffusion coefficient is given by Eq. I|66|) . As 
far as the defecton-defecton and defecton-phonon scatter- 
ing are independent, the total diffusion coefficient may be 
obtained using the Matthiessen rule. This yields 



may 

m 



D- 1 = D^ 1 = 



asA 



f T 

2^/zxa + 3a<7 -f I — 
A \0 p 



The transition from the temperature regime to the 
constant value l|39(l takes place at the temperature 



t* = r k r 1/9 



(70) 



where T k is given by the expression l|67(l . In that way the 
temperature interval may be divided into three parts. 
The region A < T < T* corresponds to impuriton- 
impuriton scattering while scattering by phonons prevails 
at T* < T < T fc . At T > T k (but below the melt- 
ing point) the classical activation diffusion comes into ef- 
fect. The diffusion coefficient behaviour is schematically 
shown in Fig. [5] The diffusion coefficient first decreases 
exponentially with decreasing temperature to Tfc, then 
increases rapidly (~ T~ 9 ), and when T < T* goes to a 
plateau which depends on the concentration (to be more 
precise, on the mean free path I ~ l/{xa)). If xa — * 1, 
then the mean free path / — > 1, and T* — ► T k . It is 
seen, therefore, that the observation of the temperature 
dependence of the coherent quantum diffusion is possible 
as far as the coherency condition Z ^> 1 is fulfilled. 

Let us consider for definiteness the case of h.c.p. 4 He 
crystal with molar volume V m w 21 cm 3 . Then the values 
of the quantities included in the expression l|69l) which 
can be measured independently are as follows 

= 26 K, p = 10/3K, a = 3, 27A, s = 320 m/s. (71) 

The quantity cto calculated with an account of the mass 
difference of helium isotopes is equal to ^7| a o — ((^4 — 
m3)/m3) 2 = 1/9. We shall limit ourselves with this value 
because an account of the change of the elastic moduli 
|37| would be based on model considerations and would 
finally give a correction small compared to the experi- 
mental error. 

Substituting (|7T)l into (|rT9")) yields the following expres- 
sion for the dependence of the diffusion coefficient on the 
temperature and concentration. 



D- 1 = 5.6 x 10 4 



xa 0.1, 
A" + A2 



(72) 



The expressions for characteristic temperatures then take 
the form: 

Tfc=4.2A 1 / 9 , T* = T k r 1 / 9 ^T k (xa)^ 9 (73) 
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TABLE I: Quantities B, G 1741 and a as functions of the 
molar volume V m and concentration x. 



X 


V m 


B.10 11 


G.10" 8 


a 


source 


% 


cm? 


cm? J s 


cm^s^K- 9 






0.05 


20.23 


0.35 


4.8 


173 


[271 


0.05 


20.42 


0.45 


7.7 


162 


f27l 


0.05 


20.62 


0.9 


4.5 


104 


f27l 


0.05 


20.84 


1.6 


10 


71 


f27l 




9u QS 


9 8 


i n 




f97l 


0.01 


20.95 


2.32 






f39l 


0.01 


21 


2.8 


10 




f27l 


0.02 


21 


2.8 


10 




f27l 


0.0242 


20.95 


2.35 






[39] 


0.0499 


20.95 


1.9 






[39] 


0.05 


21 




24 




[40] 


0.006 


20.95 


1.0 


24 




[41j 


0.12 


21 




24 




[40] 


0.25 


21 


0.81 






[14] 


0.01 


20.70 


1.1 






[39J 


0.242 


20.70 


1.2 






[22] 


0.0499 


20.70 


0.9 






f39] 



The second term in (|72[1 corresponding to the diffusion 
law H66J) allows to determine the band width A from the 
temperature dependence only . Having the value of A one 
can calculate the cross-section a from the first addend. 

In order to study the dependence on the molar volume 
V m let us write the expression (|69|) in the form: 

D^ 1 = B~ x x + G~ 4 T 9 (74) 

where 

a 2 

B = D Q x, G = 2AasA 2 dl = 4.8— A 2 ^. (75) 

y h F 

The values of the coefficients B and G at different molar 
volumes and concentrations are given in Table 1. As far 
as the dependence of 9 p on V m is known (the Griineisen 
parameter 7 is the same as for the Debye temperature). 
Eq. (fill) enables us to find the " Griineisen parameter for 
the band width" 7a : 

7A = d(ln A)/d(ln V m ) = 4 7 - 1 + X - lG (76) 

where 

lG =d(\nG)/d(\nV m ). 

Substituting here 7 = 2, 6 and calculating 7^ w 21 from 
the experimental data given in Table 1 yield 7a ~ 21. 

This value is close to that for the Griineisen parameter 
for the exchange integral in 3 He (jj « 18). Unfortu- 
nately, the experiment is not precise enough to lead to 
unique conclusions. 

Let us turn now to Eq. (|40|l and write the Griineisen 
parameter in the form 

T = dQnD)/d(]nV m ) = 2 7A -7,o + ^ « 42^- 7 „ . (77) 



The experimental values of T differs significantly even 
within the publications of one and the same group. In 
[iH ] the value 

r = 47±ll, (78) 

has been obtained while Allen, Richards and Schratter 
reported 

T = 57±12. (79) 

A conclusion may be made from a comparison of (70) 
and (71) that the interaction between impurities (char- 
acterized by the quantity Vo) depends on the molar vol- 
ume relatively small: j v ps 5. Such a dependence can 
be obtained assuming that Vq ~ ifra 3 , where Kt is the 
isothermal modulus of elasticity [4£j ■ On the other hand 
a comparison of (70) and (72) leads to a strong depen- 
dence: 7„ ~ 16. 

Let us show now the order of calculation of defecton 
characteristics on the bases of the data in Table 1. Es- 
timation will be made for V m — 21 cm 3 , x — 0.05%, 
B = 2.8 x 10~ n [23. Comparing the values for G with 
Eq. yields the band width A = 10~ 4 K. Substitut- 
ing the result obtained into the expression for B (the 
first term in (|72|l ) at given concentration one obtains 
the cross-section a = 10 4 A/(5,6 J B) = 64a 2 . Then V 
may be obtained from $Z$: V ~ aA ~ 0,65 x 10~ 2 K 
(Eq. (|21(l gives slightly larger value of the same order.) 
Making use of (|35f) one may calculate the mean free path 
/ = aj{\plxa) w 23 a and verify the coherency condition 
I ^ a. Finally, temperatures Tk and T* are calculated 
and compared with the experimental curve D(T). This 
procedure yields: 

A k 10~ 4 K, a ~ 64 a 2 , V ~ 10~ 2 K 



T k = 1.2 K, T* = 0.7(10V) 1/9 K. (80) 

If the value G = 2.4 x 10 7 reported in Refs. HJ EH is 
used, then A ps 1.5 x 10~ 4 K. Both values are within the 
experimental error. 

In regard to the concentration regions one obtains that 
the gas approximation is valid at concentrations 

x < a- 3/2 re 10~ 3 , (81) 

and the region of interacting quasiparticles is given by 
the inequalities 

10~ 3 < x < IQ- 2 . (82) 

A comparison of theoretical and experimental results 
HI M, E3 is shown in Figs. i| ??J1|7| Fig. H 
is an experimental corroboration of the ninth degree in 
the diffusion law (|66|l . Fig. 0] shows the concentration 
dependence of the diffusion coefficient in the region of the 
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low-temperature plateau. The temperature dependencies 
at different concentrations and molar volumes are shown 
in Fig. ?? and Figs. El 13 

In order to avoid any misunderstanding it is worth 
nothing that an interpolation formula analogous to l|69|) 
has been used also in 0, EHHHil : 



D~ 



n 



x+ a e 



(83) 



This expression contains two essential differences com- 
pared to (|69[) . The first one is the presence of the Debye 
temperature 9 m 88 p instead of 9 P . The second one is 
that (|83|) does not contain the cross-section a which role 
has alread y b een considered. The same error is present 
also in |2fil44| . That is why the values of the band width 
obtained in these works turned out to be smaller by more 
than 2 orders of magnitude. 

In order to describe the temperature dependence in 
the third region in Fig. [21 one has to add the classical 
expression 



D v = —lo v exp{-$„/T} 



(84) 



where the subscript v refers to vacancies because the dif- 
fusion is dominated by the vacancy-impurity exchange. 
Typical values of the exchange frequency is u> v ~ 10 s 4- 
10 9 s _1 ~ 10 -2 K. The final form of the diffusion coeffi- 
cient is therefore: 




FIG. 1: Defecton band shape in the deformation field. The 
mean free path I is determined by the scattering on the defec- 
ton band boundaries. 

D(T) cm 2 /s 




T K 



FIG. 2: Schematic representation of the diffusion coefficient 

for X3 < X2 < Xl. 



D, 



D + D„ 



2 / rp \ 9 ' 



where D is given by eq. (|69|l . and 

?Ufz 




D° v e-^/ T (85) 



ksa 2 



Co 



2.5 ±0.5) x 10 4 K.s/cm 2 



(440 ± 15) cm _2 s, 



(86) 



(87) 



be improved by further variation of the parameters, but 
this seems not reasonable having in mind the approxi- 
mations made. It is seen that the activation energy $ 
is in good agreement with what has been observed ex- 
perimentally |27| and increases with decreasing molar 
volume as it should. The band width also manifests a 
strong molar dependence changing almost two times from 
V m = 21 cm 3 to V m = 20.60 cm 3 . The corresponding 
Griineisen parameter is approximately 25 in good agree- 
ment with eq. (|76|l . 



VII. PHONON-CONTROLLED VERSUS 
VACANCY-CONTROLLED DIFFUSION 



D° = -a 2 uj v = (3 4 10) x 10~ 7 cm 2 /s 

with kg ~ the Boltzmann constant. The quantities 6 p and 
A depend on the molar volume. The Griineisen param- 
eter for 6 p is the same as for 6n. According to Edwards 
and Pandorf for 4 He & D = 7.12 x 10 4 /V^ 60 K. For 
8 V one has 



Although the theory of quantum diffusion has been 
well confirmed by the experiment some attempts to re- 
generate the idea of thermal vacancy-assisted diffusion 
have appeared recently |46l l47l El, 143. l50 | . This mech- 
anism was proposed by Locke j^j). The diffusion coeffi- 
cient in case of vacancy-impurity scattering is of the form 
JTTJ> with A = J$ A /uJv ■ 



= 0.91 x 10 4 VC 2 - 60 K. 



(89) 



We plot in Fig. [SI Fig. and Fig. [7| the dependencies 
(|f]9"|l and (|S3|l . The quantity 6 P is calculated using (|S9"|l . 
The experimental points are from Ref. 1271 The fit may 



D 



k B J- 



IV 



4 hxi v x v oiv 



(90) 



where u> v is the 4 He- vacancy exchange rate, x v is the 
vacancy concentration, and aiv is the impurity-vacancy 
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IgD- 1 




FIG. 3: The dependence of the diffusion coefficient on 
the temperature. The straight line correspond to the law 
IgD" 1 = 91gT + 6.613. 

D(x) cm 2 /s 




10- 4 icr 3 io- 2 

FIG. 4: Concentration dependence of the diffusion coeffi- 
cient at Vm = 20, 95 cm 3 in the region of the low-temperature 
plateau after |4ll| . 

cross-section (in units a 2 ). Since the vacancy concen- 
tration increases with increasing temperature as x v ~ 
e~*/ T , the diffusion coefficient varies exponentially and 
increases with decreasing temperature as 



D(T) 




l.io- 8 [ i ; 

0.5 0.75 1 1.25 1.5 1.75 2 2.25 
T -l 

FIG. 5: Temperature dependence of the diffusion coefficient 
at fixed concentration x = 5.10 -4 and a molar volumes V m = 
20.98 cm 3 . Experimental points from Ref. 1271 Solid line: 
theoretical curve for B v = 3.33 K, A = 0.9 x 10" 4 K, D = 
5.5 x 10" 8 cm 2 /s, lov = 3 x 10 9 s _1 |g3, $ = 7.1 K; dashed 
line corresponds to Eq. l9*8t . $ = 11.5 K, u v = 1.2 x 10 9 s _1 . 




0.5 0.75 1 1.25 1.5 1.75 2 2.25 
T -l 



FIG. 6: Temperature dependence of the diffusion coefficient 
at fixed concentration x = 5.10 -4 and a molar volumes V m = 
20.84 cm 3 . Experimental points from Ref. 1271 Solid line: 
theoretical curve for 6 P = 3.39 K, A = 0.9 x 10~ 4 K, 
$ = 8K; dashed line - Eq. 19811 with the same parameters as 
in Fig.EJbut $ = 11 K. 



D IV =D° Iv e*l T (91) 

with 

D% = 1L^>L = 1.77 x io-i4 cm 2 /s (92) 

where the following values of parameters were used in 
Ref.IH 

$ = 12 K, a — 3.67 A, uj v = 0.02 K, 
J 34 = 8 x 10" 7 K, a IV = 0.25 (93) 

The fit was found not satisfactory (see Ref. I2H l5ll] 1 . 

Kisvarsanyi and Sullivan (KS) [48| used the same for- 
mula with higher values of J34 

J 34 /27T = 2.3 x 10 5 Hz, i.e. J 34 = 1.1 x 10" 5 K (94) 



D(T) 



2.10" 8 




0.5 0.75 1 1.25 1.5 1.75 2 2.25 
T -l 

FIG. 7: Temperature dependence of the diffusion coefficient 
at fixed concentration x — 5.10~ 4 and a molar volumes V m = 
20.62 cm 3 . Experimental points from Ref. I27I Solid line: 
theoretical curve with $ = 9 K, D = 1.8 X 10" 8 fzft, 
6 P = 3.47 K A = 0.6 x 10 -4 K. 
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D(T) 
2.10~ 8 

1.10- 8 

7.10~ 9 
5.1CT 9 

3.10- 9 




V m = 20.62 cm 3 
x = 5 x 1CT 4 



. 5 



, 75 



1.25 1.5 1.75 
T -l 



2 2.25 



FIG. 8: The diffusion coefficient (gg) 
and different values of ui: solid line: ui 



for fixed $ = 10.5 K 
= 1.2 x 10 9 



sees that one of the authors cites a private communi- 
cation to herself ^(|. This important point, therefore, 
needs elucidation. 

The second point concerns neglecting of the impuriton- 
phonon scattering leading to Dt = G~ 1 T~ 9 (cf. Eqs. 
(1(561 . 174|) . J73J). KS have compared the reciprocal values 
of the diffusion coefficients given by expressions (|66|) and 
(|91fl and have come to the conclusion that the impuriton- 
phonon scattering is weaker by more than two orders. 
As a result, t hey have neglected D T . However, as it 
was shown in |53j ]. they have made a plain arithmetical 
mistake when evaluating the factor G. Indeed, they have 
written the factor G in the form 

G « s = ^^ 2 r m —®n = 4-14 x 10- 6 ^6 9 D 



UJD 



UJD 



dash line: 



2.4 x 10 s ; short dash line: 



s 

4.8x10" s 



long 
9 „-i 



and slightly different values of the other parameters, 
namely: 



and have used the following set of parameters 

9 D = 30K, u D = k B e D /h= 3.9 x 10 12 s _ \ 

J 34 = 2vr x 2.3 x 10 5 = 1.44 x 10 6 s" 1 



(99) 



a IV = 1.40, uj v = 1.2 x 10 9 s" 1 = 0.93 x 10~ 2 K (95) Substituting into (EJ yields 



With a = 3.27 x 10~ 8 cm this yields: 

Df v s = 1.03 x 10- 12 e*/ T cm 2 /s (96) 



Hence, one has for the diffusion coefficient 



(97) 



instead of (|69|l . The diffusion coefficient for the whole 
temperature region l|85jl is then replaced by 



D 



KS 



-1 




(Dfy 3 )- 1 )- 1 



D„ 



(98) 



with D v from IpH) . 

In their works Kisvarsanyi and Sullivan (see also 
Refs.|4ilH3) argue that eq. ffify describes well the exper- 
imental data of Allen, Richards and Schratter |27j and 
that "the overall fit for only one adjustable parameter 
is very good". The value of this parameter for V rn = 
20.62 cm 3 and x = 5 x 10~ 4 was found to be $ = 10.5 K. 
We check this in Fig. |S| (solid line). For comparison we 
show the function (|98fl for different values of u>. Un- 
fortunately, we do not see any satisfactory agreement. 
Farther on, it is confusing that in KS theory the vacancy 
activation energy $ decreases with decreasing molar vol- 
ume ($ = 11.5 K at V m = 21 cm 3 , and $ = 10.5 K at 
V m = 20.62 cm 3 ). Some other comments on their theory 
can be found in Refs. |U EH However, there are two 
points we have to discuss here. One is concerning the 
impurity- vacancy cross-section c/y = 1.40 a 2 . There 
are no details about its calculation. The only available 
information is that "the standard central force approach 
as given by Goldstein |54(" and Lennard- Jones 6-12 po- 
tential have been used. Following the reference link one 



Gj/g = 0.45 x 10" 7 cm 2 s- 1 K 9 , G KS = 2.2 x 10 7 (100) 

in perfect agreement with the values used in our works 
and confirmed by the experiment. KS have calculated 
G^-g = 6.0 x 10 -5 , i.e. more than 10 3 times larger. As 
a result, they have neglected the larger term, not the 
smaller. 

In their Reply to the Comment [5j| , KS have made 
an attempt to justify themselves by an inappropriate 
variation of the set of parameters. They changed Qd 
from 30 K to 33 K (this Debye temperature corresponds 
to V m = 19.18 cm 3 far away from the reg ion considered; 
in fact &d = 26 K for V m = 21 cm 3 45]), and J34 from 
1.5 x 10 6 s" 1 to 15.7 x 10 6 s" 1 = 1.2 x 10" 4 K (a value 
more than 10 times larger than the experimental one). 
But even with these unrealistic values they were not able 
to reach 6.0 x 10~ 5 K. 



VIII. CONCLUSION 

We presented in this work the main ideas and results on 
quantum diffusion of defects in solid helium. We showed 
that the strong temperature dependence, predicted in 
[HITH can be explained very well by the phonon-defecton 
scattering. The interpolation formula (|85fl is in the excel- 
lent agreement with the experimental data available. The 
theory allows evaluation of some important quasiparticle 
characteristic by direct comparison with the experiment. 
The exchange integral and defecton band width were es- 
timated. The fundamental role of the narrowness of en- 
ergy band width in the defecton-defecton scattering and 
the resulting giant cros-sections predicted in [Hll7| were 
confirmed. The dependence of the kinetic characteristics 



14 



on the molar volume was analyzed as well. A critical 
analysis of some other theoretical models and calcula- 
tions was made. In particular, the inconsistency of the 
vacancy-assisted diffusion at low temperatures was 
shown. 
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